Colloid-colloid and colloid-wall interactions in driven suspensions 
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We investigate the non-equilibrium fluid structure mediated forces between two colloids driven through a 
suspension of mutually non-interacting Brownian particles as well as between a colloid and a wall in station- 
ary situations. We solve the Smoluchowski equation in bispherical coordinates as well as with a method of 
reflections, both in linear approximation for small velocities and numerically for intermediate velocities, and we 
compare the results to a superposition approximation considered previously. In particular we find an enhance- 
ment of the friction (compared to the friction on an isolated particle) for two colloids driven side by side as well 
as for a colloid traveling along a wall. The friction on tailgating colloids is reduced. Colloids traveling side by 
side experience a solute induced repulsion while tailgating colloids are attracted to each other. 

PACS numbers: 61.20.Gy, 66.10.Cb, 82.70.Dd, 83.80.Rs 



I. INTRODUCTION 

Because colloidal suspensions are ubiquitous in biological 
as well as in technological systems, their non-trivial rheologi- 
cal behavior has been subject to research for a whole century 
U Tightly connected to these rheological properties is 
the internal dynamics of these suspensions, e.g., sedimenta- 
tion istl and diffusion 10. HI S 01 ■ In the context of microflu- 
idics, the interest in wall effects in confined geometries has 
increased. A lot of theoretical effort has been focused on hy- 
drodynamic interactions between solute particles and walls as 
well as among solute particles Also the collective 

dynamics of hydrodynamically interacting particles has been 
studied UHO- But only recently, the direct particle-particle 
interactions and the resulting fluid structure mediated forces in 
non-equilibrium have become subject of theoretical research 
Q. 

An intuitive picture for this so-called depletion or solvation 
force in a non-additive suspension of hard spherical Brownian 
particles has been presented in ifTill . Assume two colloidal 
particles immersed in a suspension of mutually noninteract- 
ing spherical particles, but with a hard core repulsion with the 
colloids as well as with the container walls. This leads to a 
forbidden zone around the colloids which cannot be entered 
by the Brownian particles. Once the forbidden zones of the 
two colloids (or of a colloid and a container wall) overlap, the 
osmotic pressure on the colloids becomes non-uniform, there- 
fore resulting in a net force. This force has been measured in 
various systems, see e.g., lITsi [I6IIT7I1 . 

In a driven system, the re-distribution of the spherical 
Brownian particles modifies the osmotic pressure leading to 
long ranged interactions (long as compared to depletion forces 
in equilibrium). Some aspects of this have been studied the- 
oretically for ideal Brownian particles in ifisll . The density 
of the Brownian particles near the colloids was calculated in a 
superposition approximation based on the density in the vicin- 
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FIG. 1: The two-colloids system (a) and the colloid- wall system (b). 
The colloids are driven with velocity v. The gray areas Fi and T2 
are forbidden for the center of the solute particles, the dashed lines 
represent their surfaces i9ri/2- 
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ity of a sin gle c olloid calculated in dynamic density functional 
theory lElMlll. 

In this paper, we calculate this force for a simple sys- 
tem of mutually non-interacting spherical Brownian parti- 
cles in which the DDFT reduces to a simple advected drift- 
diffusion equation, i.e., the Smoluchowski equation. We re- 
strict our analysis of the interaction between two colloids to 
two paradigmatic cases, namely colloids driven side by side 
and behind each other In addition, we consider the case of 
a colloid moving parallel to a solid wall (not discussed in 
ifisll ). In contrast to I13il , we solve the stationary Smolu- 
chowski equation in bispherical coordinates, analytically for 
small velocities and numerically for intermediate velocities. 
We additionally use a method of reflections as a second al- 
ternative to the superposition approximation and discuss the 
range of validity of the different approximations. 
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II. TRANSPORT EQUATIONS 



calculate the pressure from the local density and get 



We consider a stationary driven three component system 
consisting of a solvent, Brownian solute particles, e.g., poly- 
mer coils, and (a) two driven colloidal particles at distance d, 
see Fig.fTfa), or (b) one driven colloidal particle at a distance 
d' from the surface of the forbidden zone at a planar wall, 
see Fig. \llh). The solvent molecules are small compared to 
any other particle and the solvent is therefore approximated 
as homogeneous medium. We do not take into account the 
perturbation of the solvent flow field by the colloi dal p arti- 
cles, i.e., hydrodynamic interactions, (as studied in 12111 ) for 
three reasons. First, we want to focus on the simplest possible 
scenario, second there is no analytical solution available for 
the flow field near two spherical particles, and third, we want 
to be able to compare our results directly to fis!]. We there- 
fore assume a homogeneous solvent flow field in the frame 
comoving with the colloids, u = (— Wa:,0, — Uj) in the case 
(a) of two colloids and u = (— w^, 0, 0) for case (b), i.e., one 
colloid near a wall, see Figs. [Ha) and (b), respectively. The 
stationary density p{r) of the Brownian particles is then given 
by the Smoluchowski equation 



dp 

dt 



= V • j = -Z? A/9 + u • V/9 , 



(1) 



where D is the collective diffusion coefficient of the Brownian 
particles and j denotes their probability current density. The 
hard interaction with the colloids or the wall leads to no-flux 
boundary conditions for Eq. ^ at the surface of the forbidden 
zone dT, see Fig[T] The boundary condition for Eq. (|T) on dT 
is 



(fi-j)l 



ar 



0, 



(2) 



with the surface normal vector n pointing out of the forbidden 
zone. For colloids with non-overlapping forbidden zones, dV 
is a spherical surface of radius R, where R is the sum of the 
radii of the colloid and the Brownian particles. Far from the 
colloids the density should approach the unperturbed density 
Pa- 

As described in the following, we will solve Eq. ([T]i in bi- 
spherical coordinates, which are well adapted to the geome- 
try under consideration, as long as the forbidden zones do not 
overlap. We use a linear expansion for small velocities and we 
solve Eq. ([U numerically for larger velocities. As an alterna- 
tive we use the method of reflections, and com pare the results 
to the superposition approximation implied in II13I1 . where the 
density near the two colloids is approximated by the product 
of the single colloid densities p^^^ of two colloids with centers 
at Ri and R2, 



p(2)(r)«pU)(i._R,)pU)(r-R2)/po 



-,(1 



(3) 



The force on the colloid at Ri is given by the integral of the 
osmotic pressure over the surface of its forbidden zone dTi. 
Since there is no mutual interaction between the Brownian 
particles, we can use the equation of state of an ideal gas to 



dSksT p , 



(4) 



with surface element dS. 

The solution of Eq. ^ only depends on the dimension- 
less velocity u* ~ ^—jj-^ the Peclet number, and the di- 
mensionless distances d* = d/R and d'* ~ d'/R. Since 
Eqs. ([TJ and ^ are linear in p, p is proportional to po- 
The dimensionless force calculated from Eq. (|4|i is given by 
F* = F/{kBTpoR^). 



A. Bispherical coordinates 

Bispherical coordinates (/^i, 77, (p) are defined via 12 

a sin r] cos 4> 
cosh p — cos 77 ' 
a sin?7sin0 

y = — ' 

cosn p — cos 77 

a sinh p 
cosh p — cos 77 ' 



(5) 



with p € ] — cx3,oo[, 77 G [0, tt] and € [0, 27r]. The sur- 
face p = ±/^o is a sphere of radius a/|sinh/xo| centered 
at (0,0, ± a coth/io), p = defines the xy-plane. For the 
case (a) of two colloids (see Fig. [11^))- 5ri/2 is given by 
p = ztpf) and the distance between the centers in units of 
R is d* = 2 cosh/^o- In case (b) we choose the surface of 
the forbidden zone r2 of the wall (see Fig. [Ub)) to lie in the 
x7/-plane and the distance between the colloid and this surface 
is d'* ~ cosh pq. The radius of the spheres R fixes the scaling 
factor a = R sinh/io. 

In order to separate variables in the Laplacian in Eq. ([T]i we 
introduce /i(r) by setting 



p{r) = \/ cosh p — cos 77 h(v) + po ■ 



(6) 



Because ^/cosh p — cos r\ ^ far from the colloids (i.e., 
/i ^ 0, 77 ^ 0), the boundary condition p ^ po is fulfilled as 
long as /i(r) does not diverge for |r| 00. From Eq. ^ we 
get for h{r) 



Dk^ 



dp^ 



1 d 



sm77 — 
sm T] OT] \ orj 



dh 



1 d^h 1 



sm77^ 



-a Ux cos (j) 



■ 1 ■ 9h dh 
smnp sin77 — — h (cosh/i cos 77 — 1) -— 
op orj 



— — cosh p sin rj h - 



tan ( 



k dh 



sin 77 d(f> 



(1 — cosh/i cos 77) — — sinh/i cos 77/1 



. dh 
— smh p sm 77 -— 
dfj 



= 0, 



(7) 
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with k = cosh fi — cos 77. The first term in square brackets, 
the Laplacian from E g. dB , is diagonaHzed by the spherical 
harmonics Y™{t], cf) U2R . With n = — on dV for > 0, 
the boundary condition (|2]l becomes 



D 



— k snin un + k — — 
2 Ofi 



(— siiih fj, sin 77 cos (p) [ h + —j= 



(1 — cosh fi cos 77) ( + 



(8) 
0. 



For the two colloids in Fig. [Ha), we will only discuss the two 
cases of u being perpendicular and parallel to the axis con- 
necting the two colloids, while for the case of the colloid in 
front of the wall in Fig.[Tlb), u has to be parallel to the wall 
in order to allow for stationary solutions. 



1. Numerical solution 

The spherical harmonics Y{^{T],(f>) are eigenfunctions of 
the angular part of the Laplace equation for h{r). We there- 
fore expand h in these functions and project Eq. ^ onto 
them. Taking into account the symmetry about = 0, 
we get ^ N'^ + I iV + 1 linear ordinary differential equa- 
tions for the expansion coefficients if the expan- 
sion is truncated at I = N. The same projection proce- 



dure yields the boundary conditions for the at /i = 

zbyUQ. This set of equations is solved numerically using 
AUTO 2000 (see http://indy.cs.concordia.ca/auto/), a software 
which solves boundary value problems for systems of explicit 
ordinary differential equations by homogeneous continuation 
starting from a known solution for a special set of parameters. 
Here we started from u* = 0, for which p{r) ~ po, using 
u* as a continuation parameter. Since fc^ multiplies the first 
term in Eq. (|7J, the second derivative of each A™' (/i) appears 
in more than one equation. This set of equations has therefore 
to be solved analytically for these second derivatives, which 
restricts us to a maximum order of = 10. The case of two 
colloids driven behind each other has in addition azimuthal 
symmetry allowing us to go to = 15. Dividing Eq. (|7]i by 
before projecting onto spherical harmonics, which yields 
an explicit system of equations directly, would lead to singu- 
larities at 77 ^ and to additional numerical problems. The 
expansion in spherical harmonics converges badly for d* ^ 2 
(i.e., when the borders of the two forbidden zones get close) 
and for large u*. In the first case the interval of p for points 
outside the spheres, where the differential equation is solved, 
i.e., — /io < /i < Po, goes to zero. For large u* sharp vari- 
ations of p{r) develop near the surfaces of the colloids. For 
this reason we use this method only for 7i* < 2 and d* > 2.5. 
As an example. Fig. |2] shows contour plots of the numerical 
solutions for u* — 1. 




FIG. 2: Contour plots of the stationary density p/ po of solute parti- 
cles near two driven colloids for u* = 1 and d* — 2.51. Bright areas 
correspond to high densities, dark areas to low densities. In (a) the 
colloids are driven from left to right, the density is enhanced between 
the colloids which leads to a repelling force. In (b) they are driven 
from bottom to top, the bow wave effect in front of the rear colloid 
and the depletion behind the one in front are reduced, so the friction 
forces on the two colloids are reduced as well. 



2. Expansion in u* 

We also seek a solution for small velocities and hence write 
h{r) as a Taylor series expansion in u* , 



h{v) = hoiv) + u* hi{v) + 0{u*^) . 



(9) 



For u* ^ we require p{t) po and hence ho ~ 0. To 
first order in u*, Eq. ^ reduces to the first term, the Laplace 
equation for hi, 



d'^hi 



1 d 



dp^ ' sin 77 977 \^ ^ drj 



dhi 



1 d'^hi 1 



sm 77^ 



= 0. 
(10) 

The boundary conditions on dV depend on the direction of the 
flow field. For colloids traveling side by side we have w* = 
leading to 



— fc smli p hi + k — — 
2 op 



apo 



sinh p sin 77 cos < 



Vk 



(11) 



while for the case that one colloid follows the other we have 
?/* = and we get 



1 , ■ 1 , ,2 

— k smn p hi + k — — 

2 op 



apo ■ 



1 — cosh p cos 77 
7^ 



(12) 
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For the first case, the solution is symmetric with respect to 
/i = and in the second case antisymmetric. In both cases the 
solution is symmetric with respect to = 0. 

The solution of these equations can be expanded in spheri- 
cal harmonics, see 1I22I1 . leading to an infinite system of linear 
algebraic equations for the expansion coefficients, which is 
solved numerically after truncation at order N. As discussed 
in the previous subsection, this expansion converges more 
slowly for smaller /io, i-S-, for d* 2, so we use N — 200 
and keep d* > 2.01. 



B. Method of reflections 

This method has been used successfully to approximate hy- 
drodynamic interactions |H, I23I1 . Here we use it to approxi- 
mate the density field p(r) for the two colloids with centers 
at Ri and R2. Starting point is the solution p*^^^ for a sin- 
gle colloid, which we calculate by expanding p^^^ in spheri- 
cal harmonics. After truncating the expansion at order N we 
solve the resulting system of coupled linear ordinary differen- 
tial equations numerically using AUTO 2000. For small u* 
we get in linear order the simple analytical result 



u cos 

2 



Oiu* 



(13) 



We then immerse the second colloid into the density field 
p^^^ (r — Ri) arising from the first one. The first order density 
pi{r — R2) reads then 



/9i(r - R2) = p^^Hr - Ri) + Pai(r - R2) , 



(14) 



where Pa]^(r — R2) solves Eq. ([T]i and is chosen such that 
pi (r — R2) satisfies the boundary condition ©at |r — R2I = 
R, i.e., on the surface of the forbidden zone of colloid two. 
Re-immersing the first colloid at Ri into the density field 
pi(r — R2) yields the second order density 



P2(r - Ri) = pi(r - R2) + Pa2(r - R-l) : 



(15) 



where Pa2(r ^ R-i) adjusts the boundary condition at |r — 
Ri| = R. This procedure can be repeated and the itera- 
tion is assumed to converge to the exact solution for two col- 
loids. The approximated solutions pi,p2,--- always solve 
the Smoluchowski equation ([T]i and additionally satisfy the 
boundary condition on one of the colloids. 

In first order of u* the force on the colloids resulting from 
pi can be calculated analytically. For higher orders in reflec- 
tions as well as for large u*, the changes of coordinates from 
one colloid to the other in the reflection procedure has to be 
performed numerically. 



III. COLLOIDS DRIVEN SIDE BY SIDE 

In this section we study the case of two colloids driven side 
by side through a solution of solute particles. We start with an 
expansion in the dimensionless flow velocity u*, where all the 
introduced approaches can be evaluated analytically, and later 
discuss numerical solutions for higher velocities. 



A. Expansion in u* 

1. Bispherical coordinates 

Due to the symmetries of the system and because the right 
hand side of Eq. (fTTI) is proportional to cos (j) the expansion of 
hi simplifies to 



hi 



1=0 



C} cosh 



l + llP 



P;^ (cos 77) COS0 . (16) 



Projecting Eq. ( fTTI ) on Pj^ (cos r/) cos (j), one gets a linear sys- 
tem of equations for the . Here we need the following ex- 
pansion (see ll22ll ) 



1 



\/ cosh fi ~ cos ?/ 



= \/2 ^e~('+5)lMlp,(cos?y), (17) 



i=0 



in order to obtain with a recursion relation for the Legendre 
polynomials 



sm T] cos ( 



\J cosh p — cos ?7 



1 



1=0 



2/ + 1 



(18) 



X [-Pi-i(cosr7) - P,Vi(cos?/)] cos0, 

which is used for the projection of the right side of Eq. (fTTT i 
onto Pi (cos 77) cos (/). 

To first order in u*, the force defined in Eq. ^ in z- 
direction vanishes by symmetry, so there is no solute particle 
mediated interaction between the two colloids. A term linear 
in u* would change sign if u* changes sign. With the same 
type of argument one can show that there is no hydrodynamic 
interaction between the colloids in this setup. The friction 
force parallel to the flow u, i.e., in x-direction, is enhanced 
compared to a single colloid, since more solute particles are 
collected in front of two colloids. This force as a function 
of the distance d of the centers of the coUoids is depicted in 
Fig.S 



2. Method of reflections 

In first order in u*, Eq. ([T]i reduces to the Laplace equa- 
tion. The solutions which do not diverge for r ^ 00 can be 
expanded in terms of spherical harmonics as follows 



(19) 



The adjusted densities can be expanded in the same way. 
Starting with Eq. ( fTSl l for a single colloid, we get for the ad- 
justed density on the second one 



Pai{r - R2) = u*po 



i?2 



1 i?3 



2 4dV (r-R2)2 



cos 9 



(20) 

where 9 is the angle between v and r — R2 . The dots represent 
higher order terms in which are important for the second 



5 



2.55 
2.45 
2.35 
- 2.25 
2.15 
2.05 



□ 

\ "'a 



bispheric □ 

superposition 

reflection 1st order 

reflection 2nd order 



'■a — 



■-■B — 



2.5 



3 

d/R 



3.5 



FIG. 3: Friction force F* of ideal solute particles on one of the two 
colloids driven side by side to first order in u* . The horizontal line 
at \F*\/u* = 2jl 2.1 is the limit for d —^ oo. The same force ap- 
plies to a colloid which is driven parallel to a planar wall at distance 

d' = d/2. 
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ITG. 4: Repelling force F* of ideal solute particles between two 
colloids driven side by side with velocity u* = 1. The same force 
acts on a colloid which is driven parallel to a planar wall at distance 

d' = d/2. 



reflection but do not contribute to the friction force in first 
order The density pi on the surface of the forbidden zone of 
the second coUoid is thus given by 

pi(r-R2)||r-R2|=i? = Pa + u Po ( 2 + ) cos6' + . . . . 

(21) 

The force on the cofloid in first order of reflection is hence 
given by 

F; = -u* i?2 po ksT ^ + + 0{u*') , (22) 

which is also shown in Fig. [3] The second reflection changes 
the force to order O (^) and higher, but not the term cx 
Hence the first order reflection gives already the exact asymp- 
totic form for d ^ oo. This can be seen in Fig. [3] the graphs 
for first and second reflection have the same asymptotic be- 
havior. The second reflection shows good agreement with the 
linear bispherical results. 

B. Discussion and numerical solution 

Fig.[3]compares the drag force for small velocities obtained 
in bispherical coordinates and with the method of reflections 
with the result from the superposition approximation. With 
the solution for a single colloid in first order in u* given by 
Eq. ( fT3] l we get in the superposition approximation the drag 
force 

= ^u* PO ksT Q + ^) + 0{u*') . (23) 

By construction, both the expression obtained with the method 
of reflection in Eq. (|22] | as well as the expression above con- 
verge for d — > oo to the force on a single colloid. Both have 
the same power law dependence cx d~^, but with a different 



prefactor. For large d already the first order reflections approx- 
imate the linearized bispherical result from Sec. IIII ATI well. 
For smaller d the second order reflection result is needed. 

For higher velocities, the force between the two colloids 
is non vanishing and repulsive. Note that this is the only re- 
pulsive force (apart from possible direct interactions between 
the colloids) since there are no hydrodynamic interactions. It 
increases with decreasing distance d, see Fig.|4] The force cal- 
culated with the method of reflections to first order and the nu- 
merical result obtained in bispherical coordinates agree well 
for large d, suggesting that the first order reflections give again 
the correct asymptotic behavior for d ^ oo. 

This repelling force for fixed distance shows an interest- 
ing dependence on the velocity. For small velocities, it in- 
creases proportional to u*'^, has a maximum and goes to zero 
for u* oo, see Fig.|5] The reason is, that the force is a sym- 
metric function of u* and that the range of interaction goes to 
zero for u* oo. For large u* the thickness of the region in 
which the density is modified by the presence of the colloid 
decreases, such that the other colloid is immersed in a less 
perturbed environment. This is also the reason why first order 
reflections as well as the superposition approximation are ex- 
pected to work well. Here, for the superposition we use the 
numerical solution for a single colloid which is also used as a 
starting distribution for the method of reflections. 



IV. COLLOIDS DRIVEN BEHIND EACH OTHER 

Here we investigate the case of two colloids driven behind 
each other through the solution of Brownian particles. Again, 
we evaluate the approaches analytically in first order in u* as 
well as numerically for higher velocities. 
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FIG. 5: Repelling force F* of ideal solute particles between two 
colloids driven side by side at distance d/R — 3.1. The same force 
acts on a colloid which is driven parallel to a planar wall at distance 
d'/R = 1.54. 
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no. 6: Friction force F* induced by the solute particles on one 
of two colloids driven behind each other to first order in u* . The 
horizontal line is the limit for d — > oo, i.e., the friction force for a 
single colloid. 



A. Expansion in it* 



by 



1. Bispherical coordinates 

In the case considered here, the flow is parallel to the z- 
axis, hence u^; = in Eqs. (|7]i and (O. The system possesses 
azimuthal symmetry, so the Yi"{ri, 0) reduce to the Legendre 
polynomials Pi (cos 77). As in Eq. (|9| we write h{r) as a Taylor 
series in u* with ho = 0. hi is again the solution of Eq. ( fTOl i 
and given by 



hi = Si sinh 



1=0 



Pi{cosr]) . 



(24) 



Projecting Eq. ( fT2] i ai ji = ji^ on the P/(cos ?/), one obtains a 
linear set of equations for the coefficients 5*;, which can be 
solved after truncating at ^ = N . The term involving the 
square root of k is again projected via Eq. ( fTTI ) and a recur- 
sion relation. 

The friction force as calculated from Eq. (|4|i is reduced 
compared to the force on a single colloid, see Fig. |6] To first 
order in u*, this friction force (as well as the hydrodynamic 
drag) is identical for both colloids by symmetry. The reason 
is that p{x, y, z, u*) = p{x, y, — z, —u*) and therefore the co- 
efficient of the linear term in an expansion in u* has to be 
antisymmetric in z. 



2. Method of reflections 

The density distribution in front and behind a single colloid 
is different and therefore there exist two different first order 
solutions to the method of reflections, namely for immersing 
the second colloid in front or behind the first one. However, 
the forces on the colloids depend only on the term in the den- 
sity on the colloid proportional to cos 9. In first order of re- 
flection this term is the same for both solutions and it is given 



1 3i?3 



COS! 



The resulting force is then given by 

4 7r /I 3R^ 



f: 



-u* R poksT- 



2<P 



(25) 



©(u*^), (26) 



which is also shown in Fig.|6l Additionally one can show that 
to all orders of reflections, the force on the colloids does not 
depend on whether one starts the reflection method with the 
first or the second colloid, reflecting the fact, that in first order 
in u* , the friction forces on the two colloids are identical. As 
in the previous case of two colloids moving side by side, the 
asymptotic behavior of FJ for — > 00 is not changed by the 
second reflection. 



B. Discussion and numerical solution 

Fig.|6]compares the drag force for small velocities obtained 
in bispherical coordinates and with the method of reflections 
with the result from the superposition approximation. For the 
superposition we use again Eq. ( fT3] l keeping only terms linear 
in u* in the product of the densities. The resulting force is 
given by 

Fl = po ksT Q _ + 0{u*^) . (27) 

To first order in u*, we find the same features as in the side-by- 
side case. The drag forces obtained by the method of reflec- 
tions and the superposition approximation converge for large 
d to the force on a single colloid with the same power law 
d^^ but with a different prefactor. The first order reflection 
agrees well with the linearized bispherical results for large d, 
for small d higher order reflections are needed. 
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FIG. 7: Friction forces F* induced by solute particles on the two 
colloids driven behind each other for it* = 1. The limit of d ^ oo 
(i.e., the force on a single colloid) is given by \F* \ — 1.97. 



For higher velocities, the friction on the front colloid is 
larger than on the rear one which is shielded, see Fig. [T] For 
large d the forces obtained with first order reflection agree 
well with the forces calculated from the numerical solution in 
bispherical coordinates. Here, in order to calculate the force 
on the rear colloid, we start with the colloid in front and vice 
versa. For smaller d, the agreement is also good for the forces 
on the rear colloid, but for the front colloid higher order re- 
flections are necessary. The reason is, that the rear colloid 
influences the density distribution on the front colloid much 
less than vice versa. 

The fact that the friction force on both colloids is reduced 
as compared to the case of a single colloid can be understood 
in the following way. The bow wave of particles accumulated 
in front of the rear colloid push the front colloid, while the 
rear colloid moves in a region of reduced particle density and 
experiences less friction. 



V. COLLOID NEAR A WALL 

Here we consider the case of a colloid, which is driven par- 
allel to a planar wall. As mentioned above, the surface /^j = 
is the aiy-plane which we choose to be the surface of the for- 
bidden zone in front of the wall. So the boundary condition 
Eq. ([8]) must be evaluated at /i = /io and /i = 0. The flow 
field u is parallel to the wall, u = u^Bx and Eq. ^ at fj, = 
becomes 



dh 
9/i 



= 0. 



(28) 



p,=0 



This condition is fulfilled by the solution for the side-by-side 
case since it is symmetric about /i = 0. The condition at ^ = 
/io is also satisfied. So the solutions for these two cases are 
identical including the forces on the colloid, i.e., it is repelled 
by the wall and the friction is enhanced as compared to the 
case without wall, see Figs.[3]to|5] This is reminiscent of the 
method of images used in electrostatics. Fig.[8]shows contour 
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FIG. 8: Contour plots of the stationary density p/ po of ideal solute 
particles for a colloid driven from left to right parallel to a wall at 
distance (a) d'/R = 1.26 and (b) d'/R = 1.54. (a) is exactly the 
upper half of Fig. |2ta). 



plots of the stationary density distribution for two values of 



VL DISCUSSION 

We find repelling solute particle mediated forces between 
two colloids driven side by side which vanishes to first order 
in the velocity. This force as a function of the drift velocity 
shows a maximum and goes to zero at infinity since the range 
of the interactions goes to zero in this limit. The friction on the 
colloids is enhanced compared to a single one. Via the method 
of images the same applies to a colloid which is driven parallel 
to a planar wall, i.e., it is repelled and its friction is enhanced. 
When the colloids are driven behind each other, the friction 
is reduced and to first order in the velocity identical for both 
colloids. For higher velocities, the first colloid feels a larger 
friction than the second one. 

For the behind-each-other case, we find reasonable agree- 
ment between the results in bispherical coordinates and the 
superposition approximation. For the side-by-side case the 
agreement is less good. Our explanation is causality. Due to 
the motion of the colloids, perturbations of the density behind 
a colloid hardly influence the density around it, but perturba- 
tions of the particle density in front of the colloid do influence 
the density on its surface significantly. When the colloids are 
driven behind each other, the first one is almost decoupled 
from the second one and back-coupling effects between the 
colloids, which are neglected by the superposition approxi- 
mation, are less important. So this approximation works well 
in this case. For the same reason, the first order reflection on 
the rear colloid already agrees very well with the bispherical 
results for w* ~ 1. For the side-by-side case, back-coupling 
effects are more important since the colloids are on the same 
level and the superposition approximation gives worse results. 
The flow limited persistence time for one of the solute parti- 
cles between the two colloids is given by t = -j^. For large 
velocities this time goes to zero and a particle has no time to 
diffuse from one colloid to the other Therefore back-coupling 
effects become negligible. The consequences of this can be 
observed in Fig.|5] the superposition-results approach the first 
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order reflection in this limit. 

Comparing the forces on the two colloids driven through a 
suspension of ideal Brownian particles obtained in bispherical 
coordinates, with the method of reflections, and in the super- 
position approximation we find a qualitative agreement be- 
tween all three methods in the range in which they can be ap- 
plied. All methods are applicable for large distances, however, 
the superposition approximation yields the wrong prefactor in 
the asymptotics. For large u* the numerical algorithm using 
bispherical coordinates shows bad convergence but the super- 
position approximation and the method of reflections show 
the same qualitative behavior In summary we find that the 
method of reflections is more reliable than the superposition 
approximation. The numerical solution using bispherical co- 
ordinates has only a restricted range of applicability. In or- 



der to further test the method of reflections a finite element 
scheme could be used, which would allow for localized refine- 
ment of the discretization right at the surface of the forbidden 
zones. This would also allow to study the dynamics of the 
system. 



Acknowledgments 

We thank S. Dietrich for support and fruitful discussions. 
M.R. acknowledges funding by the Deutsche Forschungsge- 
meinschaft within the priority program SPP 1 164 "Micro- and 
Nanofluidics" under grant number RA 1061/2-1. 



[1] A. Einstein, Annalen der Physik 19, 289 (1906). 

[2] A. Einstein, Annalen der Physik 34, 1911 (1911). 

[3] R. Klein and G. Nagele, Nuovo Cimento D 16, 963 (1994). 

[4] G. S. UUmann, K. UUmann, R. M. Lindner, and G. D. J. 

Phillies, J. Phys. Chem. 89, 692 (1985). 
[5] G. D. J. Phillies, G. S. UUmann, K. UUmann, and T.-H. Lin, J. 

Chem. Phys. 82, 5242 (1985). 
[6] G. D. J. Phillies, C. Malone, K. UUmann, G. S. UUmann, 

J. Rollings, and L.-P Yu, Macromolecules 20, 2280 (1987). 
[7] R. Klein and G. Nagele, Current Opinion in Colloid & Interface 

Science 1, 4 (1996). 
[8] I. K. G. Dhont, An Introduction to Dynamics of Colloids, vol. II 

of Studies in Interface Science (Elsevier, Amsterdam, 1997). 
[9] J. Happel and H. Brenner, Low Reynolds Number Hydrodynam- 
ics (Prentice Hall, New Jersey, 1965). 
[10] S. Kim and S. J. Karrila, Microhydrodynamics : principles and 

selected applications (Butterworth-Heinemann, Boston, 1991). 
[11] S. Harris, J. Phys. A: Math. Gen. 9, 1895 (1976). 
[12] B. U. Felderhof, J. Phys. A: Math. Gen. 11, 929 (1978), times 

Cited: 269. 

[13] J. Dzubiella, H. Lowen, and C. N. Likos, Phys. Rev. Lett. 91, 



248301 (2003), cond-mat/0306069. 

[14] S. Asakura and R Oosawa, J. Chem. Phys. 22, 1255 (1954). 

[15] C. Bechinger, D. Rudhardt, P. Leiderer, R. Roth, and S. Diet- 
rich, Phys. Rev. Lett. 83, 3960 (1999). 

[16] L. Helden, R. Roth, G. H. Koenderink, P. Leiderer, and 
C. Bechinger, Phys. Rev. Lett. 90, 048301 (2003). 

[17] J. Crocker, J. Matteo, A. Dinsmore, and A. Yodh, Phys. Rev. 
Lett. 82, 4352 (1999). 

[18] U. M. B. Marconi and P Tarazona, J. Chem. Phys. 110, 8032 
(1999). 

[19] U. M. B. Marconi and P. Tarazona, J. Phys.: Condens. Matter 
12, A413 (2000). 

[20] A. J. Archer and M. Rauscher, J. Phys. A: Math. Gen. 37, 9325 
(2004). 

[21] M. Rauscher, M. Ki'iiger, A. Dominguez, and F. Penna (2006), 
in preparation. 

[22] P. M. Morse and H. Feshbach, Methods of theoretical physics 

(McGraw-Hill, New York 1953). 
[23] J. L. Aguirre and J. T. Murphy, J. Chem. Phys. 59, 1833 (1973). 



